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SHOPSIS 

V. * 

A numerical analysis of steady, normal shock wave structure 
in binary gas mixtures of inert, monatomic gases is made. The 
analysis is based on the moment method of kinetic theory coupled 
with Mott-Smith ansatz. . The governing equations of this two 
temperature- two velocity model are in no way restricted by initia 
concentration of the species, mass ratio of the species or Mach 
number. The system of equations which describes the structure 
includes the conservation of mass of each species, the conservatio 
of momentum and energy of the mixture and two transfer equations. 
The non dimensional parameters which appear are the mass ratio, 
free stream Mach number and initial concentration ratio. 

Only the .rigid sphere gas and the Maxwell molecule have been 
examined to study the effect of molecular model on the structure. 
Regarding the transfer equation, the choice had been dictated by 
simplicity ard is either or moment of Boltzmann equation 

where is the velocity of the particle in the x-direction. Bor 
particular,: cases, a given type of structure suggests more than 
one scale for a complete transition. Gomoarison with existing 
experimental and theoretical results is made. 

The governing equations form a system of two coupled, 
nonlinear, autonomous, ordinary differential equations whose 
boundary and regularity conditions are prescribed at the two 
infinities. This has been reduced to an initial value problem 
by analysing the singular points in the phase plane of the system 
and the existence and uniqueness of the solution have been 
considered. 



Finally, the axhitxariness of mament method has been 
discussed and a criterion to remove the arbitrariness has been 
suggested. For a special case of Mach number tending to infinity 
and rigid sphere gas numerical result has been presented which is 
nearest to the true Boltzmann solution in Ig norm when the 
assumed representation of the distribution function is bimodal . 




INTRODUCTION. 


The problem of the structure of a shock wave is 
one of the classical problems in the sense that it is 
conceptually simple, physically inp~rtant and has no 
boundaries - boundary conditions are put at _+ oo ; still 
it is a challenging problem. 

The equations of con^ervf’tion of Eiass, momentum 
and energy for a one dimensional steady flow of an 
ideal gas admit a non-trivial solution which represents 
a possible jump in the flow quantities like velocity, 
temperature, pressure, entropy etc provided the upstream 
Mach number is greater than one. In this hydrodynamic 
description, shock wave is a plane of discontinuity. 

Jump in entropy indicates a dissipation of mechanical 
energy and an irreversible conversion of mechanical 
energy into heat. 

The consideration of steady, one dimensional flow 
of a viscous, heat conducting fluid, approaching finite 
limit values at x = - oo and x = + oo (these limiting 
values being interrelated) includes the two dissi-pating 
mechanisms viz, viscosity and heat conduction attribut- 
able to the molecular structure of the fluid. These 
mechanisms account for the irreversible transfer of 
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mechanical energy into heat and replace the plane of 
discontinuity in the flow variable by a thin layer, 
sometimes called the shock layer in which steep gradi- 
ents appear. It is believed that the two dissipative 
mechanisms play entirely different roles. Viscous 
mechanism causes the scattering of the directed mo- 
mentum of the incident gas and the conversion of the 
kinetic energy of the directed molecular notion into 
the kinetic energy of the random motion, Role of the 
heat conduction is secondary in as much as it indirectly 
affects the conversion of the mechanical energy as a re- 
sult of redistribution of pressure. (Sel'dovich and 
Raizer (1967)). 

The solution of the one dimensional hydrodynamic 
equations neglecting viscosity, but retaining heat con- 
duction was ■■■§ studied by Raylagh (l910)who proved 
that a continuous solution for the flow variable is 
possible only for weak shocks. On the other hand, con- 
itinuous solution of all flow parameters for all shock 
strengths e^dsts hy considering viscosity even if heat 
conduction is neglected, Qualitatively, the problem is 
studied hy Raylegh, IS^eyl (1949) and Gilbarg (1951). For 
a fairly general set of perfect gases for which a relat- 
ion among coefficient of viscosity, thermal conductivity 
and specific heat at constant pressure is postulated, 



explicit integrations of the equations of motion has 
been carried out, (Becker (1922), Puckett and Stewart 
(1950 ), Morduchow and Lihhy (1949)), 


As the characteristic length scale of the shock 
layer is the mean free path, the question of appropriate- 
ness of continuum and kinetic description arises naturally 
The Navier-Stokes equations provide a solution of the 
shock structure for any strength of the shock where as 
the Burnett equations fail to give a shock structure for 
an upstream Mach number above 2,2 as found by Zoller 
(l95l). This, along with the experimental work of Sherman 
and Talbot (i960), leg to a wide spread use of the Navier- 

Stokes theory for the study of the shock wave structure, 
-out 


wit^ elaborating restrictions on the Mach number. However 
except for low Mach number, the agreement between the 
Navier-Stokes results and experimental results is not very 
satisfactory. Fence, though one may admit {wm 
IpHHHBl) adequacy of the Navier-Stokes theory for low Mach 
numbers, in general, it is untenable to think of any con- 
tinuum theory as apnrooriate for a situation in which the 
mean free path is a significant fraction of the shock 
thickness, the only macroscopic length relevant to the 
problem. Consequently attempts have been made to examine 
the structure of the shock wave from the kinetic theory 
Doint of view using the Boltzmann equation'?. 
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Tbe study of the shock structure my, at first 
sight, apoear to he of academic interest only in as 
much as the downstreara conditions can he obtained 

without determining its structure and its thickness 

is very small. In fact, to the order of quantities 
studied in the classic’ll shock structure, its knowledge 
throws some light on the physical nature of the shock 
wave. One field wheri the knowledge of shock structure 
may have a direct application is the study of flames. 
However, when the global problem for flow past an ob- 
stacle is studied along with the higher order effects 
(for a curved shock wave), the knowledge of the shock 
structure becomes a necessity to know the modified 
Rankine Fugoniot conditions. In the case of corres- 
ponding mixture problem, the concentration ratio gets 
modified which cannot be determined without knowing 

t • 

the shock structure (Oberai (1965)), 

Boltzmann Equation; 


The governing equation for non-uniform flow of a 
single gas in terms of one particle distribution funct- 
ion is the Boltzmann equ'’’tion. 


Df _ c>f 

* at 


+ c. 



M + X . ^ = j 

dx m ajs 

- f f^)g dii d JD ^ 


(f,f) 


( 0 . 1 ) 
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where c, are the velocities of colliding pair of 
molecules before collision and e', c^ ' are their 
velocities after collision, and where 

f^ = f(t, X, JC^) 

f’ = f(t, s, c,*) 

fl » = f (t, X, c^) 

S - “rS 

s’ = -s' 

and dix = b db d^ (see fig. 1) 
provided 

(a) The consideration is restricted to binary collis- 
ions - i.e.j^dilute gas 

(b) F2(t, x^, c^, Xg, £ 2 ) = f (t, x^,c^) f(t,X2,£g) 

where F 2 (t, x^, X 2 , £ 2 ^ probability of 

siraaltaneously Observing 2 molecules in the state 

(2l» S±y ^2’ -^2^ • 

(This is knoxiFn as postulate of molecular chaos) 
and finally 

(c) The collision between the pairs of molecules of 

arbitrary impact parameters b, are equally probable j 
i.e.jthe function f (t, x^, does not change over dis- 

tances of the order of the collision cross section. 

The validity of the Boltzmann equation is considered 
by authors like Yvon (1935), Bogolyubov (1946), Born and 
Green (l946), Kirkwood (1947) and Grad (1958), 
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From ftie definition of f , the number density n is 
equ^l to 

n(t, x) = Jf(t, X, do (0.2) 

where the integration is perfomed over all molecular 
velocities, Similarly, the mean velocity of the mole- 
cules, the stress tensor anc’ the energy flux are defined 
by the relations 


u(t, x) = 

1 

n. 

f c f (t, X, 

0.) ^0 

(0.3) 

Pij 

h 

f (t, X, 

c) dc 

(0.4) 

ra 

^i = j- 

I 

f(t, 

X, jc) d^ 

(0.5) 


Here C = c - u is the ttermal or peculiar velocity of a 
molecule. 

Temperature is identified as the mean energy of 
thermal motion by the relation 

I kT = i f (t, X, dg. (0.6) 

where k is the Boltzmann constant. 

Methods in Kinetic Theory to solve the problem of 
Shock Structure . 

Let the x - axis be directed along the flow. Let 
^ ^ ^ 

, T^and U^j, be the mean velocity, density and 

temperature of the gas, respectively in front of and be- 
hind the shock wave. 


The gas is In equilibrium at Infinity in both 
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directions (i.e. before md after the shock), therefore, 


the corresponding distribution functions are Maxwellian; 

w3/2 


f(x-^- oo) = f^( 2 .) = N”(27^RTviexp 

-3/2 


(s. “ i 

2RTv^ 
(0.7a)' 


exp 


f(x-#>f 00 ) = fgCa) 


N4(27[RT^) 


2RTj 




(0.7b) 


The investigation of the structure of the shock wave 

requires the solution of equation (O.l) with the boundary 

conditions 0.7(a) and 0.7(b) where given 

^ ■¥: 

determination of N^, and is a p?. rt of the problem. 

There is as yet no exact solution for this problem. 
Monte Carlo trials have been made by Bird (1968) in which 
a mathematical simulation of the phenomenon is carried 
out on high speed computer. The molecular model and the 
collision dynamics are assumed and the motion of one or 
several of the chosen particles is traced by the computer. 
Such methods, whatever their merit, are limited by the 
storage capacity of the computer. 


Chorin (1972) has presented a numerical solution of 
the problem recently. The method takes formidable time 
on the computer, needs large storage capacity and places 
particularly heavy demands on the accuracy of the 
algorithm used. 
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Model equations instead of the Boltzmann equation 
have also been solved to throw light on the problem. 

In these procedures, the collision t’erra of the Boltzmann 
equation is replaeed by suitable simpler expression such 
that the gross features of collision terms are retained. 

A well linoiTO model is the Xrook’s model (1954) and the 
corresponding nroblem has been solved numerically by 
Liepnann, Narasimha and Ghahine (1962). The success of 
their method is not restricted by the value of the up- 
stream Mach number. It seems that for strong shocks 
the upstream side is ouch sraoother than is predicted 
by Navier-Stokes theory, but it is not known exactly 
whether it is a property of the Krook equation or happens 
as well for the full Boltzmann equation. 

First satisfactory, approximate solution in analytic 
form was given by Mott-Smith (l95l) for all Mach numbers. 
He postulates an appropriate form for f, containing a cer- 
tain number of parameters (which may be constant or 
function!' of x). These parameters are determined by 
taking suit>^ble moments of the Boltzmann equation. The 
distribution function, according to his postulate, is of 
the following form; 

3 

“2 

f(x, ^ « N’*(x) (27TRT) exp 


(S - m.) 

2RT 
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4* N (x) (2 7TRT ) exp 


(c - F*i)^ 

/-w f 


2F.T 


( 0 * 8 ?^) 


witli the boundary conditions 


Lin 

X-^ - CO 


N“(x) = Nq, = 0 (0.8b) 


The quantities U and T are respectively the upstrean 

velocity and tenperature. The four unknowns N’’(x), 

^ ^ ^ 

N (x), U and T are deternined with the help of three 

conservation equations (viz. mass, nonentun and energy) 

and one (arbitrary) nonent equation. The final answer 

will depend upon the particular extra nonent chosen and 

the nolecular nodel assuned . Despite these linitations, 

I-Jolway (196^”?) while exanining the effect of collision 

rjodels upon the shock wave concludes that the Mott-Snith 

solution is a mich better approxination to the solution 

of the BoltziHann equation than that of the Krook kinetic 

equation. The nuia.erlcal computation by Chorin (l97:jSSj) 

also lends credence to the results of the Mott-Snith 

theory. The bimodal character of the distribution 

function is also confirned by the calculations of Bird* 

Salwen et al (l964) attempted to improve the Mott-Snith 

solution by taking multimodal distribution. It nay be 

* 

U 

mentioned here that the ratio ^ being constant is an 

u 

essential feature of Mott-Snith ansatz (for the single 
gas case) and we shall refer to it in the mixture problem. 
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■ -As far as the prohlen of shock rrave structure la 
gaseous mixture is coii«erned, motivated hy the success 
of Mott-Snith theory, its simple extension to the case 
of gaseous mixture was first given by Oberai (1965). 
(Fujimoto (1965) also attempted the solution independ- 
ently but he did not consider the conservation of mo- 
mentum and energy in the shock and consequently his re- 
sults apceared to be unrealistic.) Oberai (1965) had 
given the governing equations for shock wave in a gase- 
ous mixture but he presented solutions only for small 
concentration ratio. The governing equations are the 
following. 




•^ 3/2 

^i £.) = ”iW .(27TRjT) 


exp < 




-3/2 

n^(x) (aJTRj^T?) exp 



y 


i = 1*2 (0.9) 

t 

f j s are the distribution functions for different 
species. 


Writing nj(x) = nT(x) and 
nt(x) = N — nJ" (x) 

u y 1 


dN. (x) dNt (x) 

_ . ^ = 0. 1 = 1,2 


( 0 . 10 ) 
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mg ^ (U^ + RgT) - 


(U +R 2 T") 


*\ U ) ^^2 


7 K 


= 0 ( 0 . 11 ) 


■4 


m. 4 0(0*^ + 5R,T) - 


* , * 

U (U + 5R,T 


*\ V I 


m„ I U(U^ + 5RpT) - U*(U* i 5RgT^) \} • 

u 3 


^ ^ -‘2 


ar ■ ® 




€ (1 +€) 1 - 


3 1 + OC^ 

5 d 


(0.12) 


a-*) * ®Xj2 (- 

1 4 ^-{^2 {ff^ - “2 ( tt ^ - € )j "1 

^ ( .» 4 ^(1 +<f) 1 + OC^ 

" 1 ^ irTT? (i-p)" 


2A. 


(4^ '■■•>} 


Ng + ^2 (1 -C ) Nj + 



X ftj? 

10b 


X, 


12 


{«, * «2(T-y 


•% Nj 


(0.13) 


dN, 


3 1 t OC0 I 3 

e (1 + 4 ) j 1 - 5 -p ;i -iX) J 3^ 

4 


^^12 ) '^■^1 


1.A_ ^4.,^ 13^ ^ . 


2A, 


(1 +^) 
1 - € 


5 1 


(1 +py 


(1 + 


5 1 " 


i- 


1 + € 


€0X22 (1 -i) * Xj2 (-^1 JTTf)^ " ^ (1 


1 * (cfi 1 - ^ 
rnr " ^*2 »-. „>2. 


(1 *^y 


^2 * *22 A 2 (1 -e ) N 2 + 


X 


-4A, (1 -P) (1 -6)_ ^ ^ 


12 


(1 +p)' 


(1 +j3) 


^1 


Nj Ng 


(0,14) 


Where 6 = T ( 1 + ”2 ^ 

M 


and M 


= U^/ ( I Op T) 
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and the boundary conditions are 


Lim Nj = 1, = 


0 


= oc m; = 0 


(0.15) 


Tanji’enbaurs and Scott (1966) pointed out that the 
forfeulatlon of the above nentioned problen vas inconsis- 
tent in as much as the constancy of U and T inplied 
that component gases would not separate, Tanja'enbaun and 
Scott concluded th«t the mixture problem could not be 
solved using bi-modal description. 

In view of the above we shall discuss in chapter I 
the formulation of the following problem for the mixture; 

I, I'lultinodal distribution function 
II, Biraodal distribution function of the components 
with suitable modification. 

The discussion of the equations of problem II and 
their solution comprise chapter II which is the main 
body of the thesis. 

As has been said earlier the shock wave struetme 
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in the Mott-Sisith theory depends upon the particular 
extra moment chosen and this arbitrariness is one of 
the shortcomings of this method. The study of this ar- 
bitrariness was initiated by Mott-Smith himself and he 

2 3 

nresented two solutions correspondlag to q and C 

X X 

moments* Rosen (1954) attempted to remo’we this arbi- 
trariness by means of a restricted variational princi- 
ple, Gustafson (i960) showed that the restricted varia- 
tional equation used by Rosen was a special case of a 
transfer equation and Oberai (1967) pointed out the non- 
uniqueness of Rosen's method. Later Narasinfiir (1966 ) 
proved that to every restricted variational principle 
there corresponds a transfer equation. However, one 
way to remove the arbitrariness is to find out that 
bimodal distribution which has a minimum deviation, in 
the least square sense, from the exact solution which, 
of course, is unknown (Oberai (1967)). However, this 
procedure necessitates the calculation of the collision 
Integra’! itself, instead of its moments, Oberai pre- 
sented solution only for the special e^.se of and 

rigid sphere gas because in this case the upstream Max- 
wellian becomes a delta function and integration in the 
velocity space can be performed directly. Later, 
Narasimha and Deshpande (1969 ) obtained shock wave 
structure for the case of arbitrary Mach number using 
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this criterion by expressing the collision integral 
in closed form and evaluating the integrals needed in 
the foriaulation of the governing equations nuaerically. 

The criterion to renove the arbitrariness can be 
applied to the mixture problem. In chapter III we shall 
discuss the formulation of the mixture problem for the 
case of arbitrary Mach number and present the solution 
for the special case of oo and rigid sphere molecules. 


***** 



ClAPflB I 



In this chepter we shall formulate the problem for 
a binary gas mixture retaining the essential features of 
the Mott-Smith ansatz i simple minded approach will be 
to write the distribution functions for different gases 
as follows: 

f (c - U i)^ j 

(x) (2AEiT) exp|- + 

+ ~ I ( i)^) 

(x) (27 VRj^T*) exp< --. ^2 " r T*"^" 

. i i 


for i = 1, 2 


(1,1 a) 
b 


U and T are upstream velocity and temperature. IP- and 

T*, though unknowns are postulated to be constants. 

- + - 

So that in all, there are six unknowns: N^, N^, 

+ 

Ng, U* and T* and the boundary conditions are the 


following : 


■•oo, y* “ 1, n'q 


"l = f'2 = 0 


(l.2a) 


X— ^ + c<; “ ^2 “ ^ 


(1.2b) 


Boundary conditions (l , 2b) imply that U* and T* are 
the downstream values of velocity and temperature. 
This was done by Oberai in 1965, 


The conservation equations are the following: 


UNj + J * ® 

[ un; . . 0 


(1.3) 


(1.4) 
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a m 

t 

r 2 


dx I 


d 

dx 


1 1 h 

+ R^T) + (iff + R.T) 

fl 

O 

(1.5) 

r- 2 

1 

* 

_^U 2 IT 

+ 5RjT) + Nj^(U* +'5 RjT) 

n 

o 

(1.6) 


Integrating tfee above e(^uations (1.3) - (1.6) and 
evaluating the constants of Integration from the up~ 
stream conditions, we obtain 


= 1 


Ng = OC 


TJ* 

'M 

U 1 


2: 

V 


N, 


(1.7) 

(1.8) 


c 2 + — - 2 —7 

m.N. |u (U* + R.T*) - U*(U^ + RiT^ = 0 (1.9) 


^ “i^i L(U*+5Rj^T*) - (U^ + 5E^T)^ = 0 (l.lO) 


Equations (1.9) and (l.lO) can be solved for U* and 
T*. Ignoring the trivial solution 
U = U* and 
T = T*, 

we get 


„* . 5 R,T 

£_ = 1 (1+ I - -i ^ } 

tr 4 ^2 




and P = 1 + 


1 - 


' /U* 



<u 


K ^^^2 


( 1 . 11 ) 

( 112 ) 



XO' 


E, Kp 

Where fi = ^ - TT' > 

} Eg 

Is the molecular mass of the first gas and mg Is 
the molecular mass of the second gas. 


As U* and T* have already been ostulated to be 
constants equations (l»ll) and (1,12) imply that the 
expression ' 

^ 

. . '.-■•nJ- i-PnJ 

does not depend on x throughout the shock wave - and 
hence is a const'^nt. Equations (i,7) and (l,8) imply 
the relation 



and the boundary condition (l.2b) implies that 


Lt 

X-» CO 


N. 

N, 


(def) 


m 


1 

oc 


The constant v«lue of the expression 

K * ^2 


+ 


n; 


must equal to 


N J + N| 

N| + /5N* 


1 + OC 
1 + OC 




or 


:NJ 


= OC 


(1.13) 
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(l-13) along with (1.7) and (1*8) implies 



throughout the shock wave so that the ratio of the 

— -f 

gas densities viz, ^'^2 ^2 equal CC for all points 

n“ + 

of the shockwave. This was pointed out by 

and Scott (1966 ) who concluded th^t the mixture problem 

cannot be solved via. a bimodal ansatz. 

'^e shall examine whether the above shortcoming can 
be overcome by assuming distribution function to be 
trimodal in character. 


¥e assume 


^1 ^ ^*11 ^12 ^13 

(1.14) 

^2 " ^21 ^^'22 ^ ^'"23 

(1.15) 

where = M(N”,u , T, Rj^^) 



Mi 2 = M(Nj,li*', T*, R^) 

M^3 =K(n^, U*, T*, Eg) 

= M(n;, U, T, Rg) 

Mgg = M(N2, R*, T*,. Rg) 

Mgg = ^(ng, tr*, T*,, R^) 
where M(F, IT, T, R) denotes 
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-3/2 

N(x) (2 7?HT) exp 

The choice of and M^g is dictated by the considerat- 

ion of simplicity, l,e., the problem should remain tract- 
able. There are 8 uc’CTowns namely n^, Ng , 

+ * * 

Ng , U and T and the boundary cpnditions are the 
following, 

x^- 00 , n“ = 1, Ng = OC , = Ng = n^ = Ug = 0 

X-*. + 00 , N” = Ng =. n^ = ng = 0 


The derivation of the conservation equations is a 
standard one (see for example Oberai (1965)) and only the 
final results are given below: 




= 0 (1.18) 
U (U^+ 5R^T) + U*(U* + 5R^T*) + n^U* 

mgJ^Ng U (u^ + 5RgT) + Ng U* (U* + SRgT*) + nglJ* 

(U""^ + 5R^T*)j| = 



0 


(1.19) 
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* 


If, after integreting the system of equations 
(l,16) to (1.19) and evaluating the constants of inte- 

u' 

gration from the upstream conditions, one expresses ^ 
* 

T — — 

and in terms of and etc. one gets 


1 + OC - N, - N. 

1 £ 


* H T 

IL = i + 1 -111 

U 4 4: TT^ -I /Vi TbiT*" i* 

U 1 + OC^- 


T 


* 


( 1 . 20 ) 




(1 - N ") (1 + 0 C + (1 -^^) (1 + 

and ~ ‘1 


If (1 + 0C|S ) J 

( 1 . 21 ) 


-here ^ - i + i ^ j-..-±„Qg. 

.here t-4+4 ^2 1+ OC^ 


and ^ = m 


Equations (l,20) and (l,2l) explicitly show that for 
the formulation to he consistent, two extra conditions 
are needed jnamely| 


+ ''1 

and N„ = 

2 ^ 


( 1 - N~) ( 1 - QDp) 



(1.22) 

(1.23) 


These are artificial restrictions (predetermining sepa 
ration) and resemble the one pointed out hy Tannenbaum 
in the bigiodal description. 
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Addition of one more Maxwellian in the distribution 
function makes the problem formidably complicated, with- 
out helping remove the shortcoming, Y!e conclude that 
multimodal description does not lead to a physically 
acceptable solution of the mixture problem. 


The second v^ay to attack the problem is to modify 
the bimodal description in such a way that the above 
mentioned inconsistency does not arise. This was done 
by Oberai in 1966, ^le nostulated the distribution 
function in the following form; 


(x) 


3 

2 


(2 a R^T) 


exo 



(Q - Ui)^ 
2R^T 


} 




nJ' (x) I" 


2 7^ R^T 


3 

^2 

(x)J exp 


(a - 

2 R^T'^ (x) 


and --4 

f„ = N“(x) (2r\R T)*- exp 


(1.24) 


Ng (3£) 


2 

3 

2.’ 


(e - 


2 RnT 


2 ^ J 


2;\RpT (x) '"i-. exp 


(g-U'^(x)^) 

aRgT'^Cx) 


2- 


(1.25) 


with the conditions 

x-^- oq n“ = 1 , n“ = OC , nJ = Ng = 0 (l. 25 a) 

x-j*- +oq = 0 


The following are the conservation equations when 
the distribution function f is as given in the equation 
(1.24 - 1.25) 



d 

dx 

UN^ + ' 

= 0 


**1 


d 

jUMj + u-^nJJ 

= 0 

dx 


d 

dx 

f ^ - r 

A . 2 — El. 1 1 

1 = 1 1 L. 

.1- (u' 


(1.26) 

(1.27^ 


2 “i") 

+ nT (u"^ +R.T'^) 

1 ■ I ' 1 W 


11 j 

0 ^ (1.28) 


dx 1 1 


2 !' 

.2r . m. I 


1 ■*■ 5RjT’^)j 

= 0 (1.29) 


Integr*=ting equations 1.26 - 1.29 and evaluating 
the constants of integration from the upstream condit- 


ions, we obtain 


n 7 = 1 - fj ■*■ 

^^1 ^ u 1 


^2 ^ “ U ^2 


(1.30) 

(1.31) 


521 m,- N7ru(u'^ + R.T"^) - u'^(u^ + R.T) = 0 (1,32) 

• ji X X* f X. X 

J_ _ 1 l_ 


m.Nt (u'^ + SR.T"^) - (U^ + 5R,.T)”1 = 0 (l.33) 

^ ^ ^ ^ -J 

(1,32) (1.33) ^ ^ 

Equations A and/ce.n be solved for U and T . Ignoring 


+ '• + 

the trivial solution U = U and T = T , we get 

u" 1 _ . /: '■I " 1 

« ^ L 1 tj2 nJ . 

tI , s! r, "1 ”i 

T - ^ SRjT I I \ n;* Kj 


(1.34) 


(1.85) 



where [}, = — - , the ratio of the molecular masses 

/ 

Using the four conservation equations, four unknowns 
Nf, Ng, U"'", T^have been expressed in terms of and 

The functions and are to he determined from addi- 

tional transfer equations. No inconsistency arises be- 
cause U%.nd t"^ are already functions of x thereby allow- 
ing for separation of the species. 


Limiting value of any function of interest at 
downstreaii" infinity will be denoted by asterisk e.g. 
denotes (x) when x approaches -+ oo 


(1.34) 

Taking the limit of the equations/- and 

(1.30) (1.31) 

x^oo and making use of and ^ , we have 


(1.3 
l_ as 


35) 


- L 

U "4 




1 + 


1 + PC 
1 + 2^3 


5R^T 


) = € (say) (1.36) 


illi* _ 

T 


1 + 


1 “C- 

4C. - 1 


(1.37) 



(1.38) 

(1.39) 


Relations 1.36 to 1,39 nay be recognised as the Rankine 
Hugoniot relations if we .‘j^serve th^t the expression for 
Cp the specific heat at constant pressure of the 
mixture, evaluated upstream is 


C 


P 


1 + PC 

2 1 + OCyS 



so that the upstreaiD I'^ach nur.her K which is defined as 

(| GpT)"^ 

is related to£as follows; 

c- = ^7 ) 

2 

Oberai took monient for Maxwell molecules to deter- 

mine the shock structure and presented the solution for 

small values of the concentration ratio^. The transfer 

2 

equations with 0 moment are the following, 

X 


d I U 




r r 

'12 I 1 l_ 


r+ + + 


^12 1 I ^2 “ u"" > 


+ 


+ R, (T^^-T) 


+ ( oc- ~ N2) (uu"" - u ) + (N2 - ocN^) T-rj 


+ (4A^ - 2ij„) i 


2^ I 1 + n 


_Hh + + 


(N^ H- OD N* - N^) 


(O - V*)^l - >*2 (1 - g- nJ') N* (U - U*)^ (1.40) 

i| { W •'J [ 3^2 - T) - 


(OC- ( UU"^- U^) + (1- ~ N^) 


nJ (UU'*'- ) + (® nJ - Np — 


R. 


^ITT^ 


N+ N^) (0 - U-") 3 


2 
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- Xgg Ag (oc- 2- Ng ) Ng (U - (1,41) 



and A^, Ag are certain Integrals for Maxwell molecules 
discussed in tlie appendix, 

3 

Fe also give here below the C moment transfer 

iX! 

equations for Maxwell molecules; 


d 

d?^ 




* p op 4: 2 

N“ ( U^ + eU^Rji^T + T ) + + 


2 2 
3 T"^ 


2 

-^23 (U’‘’ + 3R.jT )| + 

j^U^(U'^ + R^T'*’) + U^(U^ + 3Rj^T) - 2UU‘^(U^+R^T)J + 


Z4 

Z5 


U^(2Rj 


T + 2R^T‘^) i + S3 


(U"^ + R.T'^) -■ 


UU'*' (U'^ + 3F.j|_T'^) 
h S4 UU"^ (U^ + R^T) + UU'^ (U"^ + SRjlT'^) - 2U^ (U"^ +R^T‘^ 




Z5 I UU'*' (2R^T + 2R^T'^ 


4- 



N“ 


+ x^g |azi 


2 ? 2 
UU'^ (U* + 3RoT'^) - 3U^ (u"^ + RoT"^) + 





3UU'" (U^ + R^T) - + 3R^T)"] ^ P.Z2 [<UU'"- U^) 

(2E2'r'^ + 2F.^T)'] + ilZ3 _UU^(U^ + + ^^1^^ 

+ AZ 4 j^ U^ (U"^ + U^(U^+ 3 R^T) - 2 U'^U (U^+R^T) 

+ AZsfu^ ( 2 R 2 T'^ + 2 F,^T)"[|it“ + 


^12 


'/Zl fu^ (U^ + 3 RoT) - 3U'^'J (u^ + RpT) + 3U^ 


L 


+ AZ2 


( 2 R 2 T + 2 R^T'^) 


(U’^ + R^T"^) - U'^U (U"^ + 3 R^T^) 

(u^ - u'*'u)J + : 

AZ 3 (U"^ + R^T"^) - U'^U (U"^ + 3 P.^T'^)'^ 

o 

AZ4 j^U'^U (U^+Rg^) ■" + ^t"") 

+ AZ5 ( 2 R^T + 2 R^T'^^| N^ (x) n” (x) (1.42) 


and 

±_ 

L 


Ng (U^ + 6 U^ RgT + SRg T^) + N^ (u'^ - 1 - SU"^ R 2 T'^ 


+ 3 R 2 T 


"12 

.2 


2 — j 

^ )> 

^Bzi pau'^ (u'^ 


+ 3 RpT'^) - 3 U^ (U"^ + RpT"^) + . 3 inj'^ 


(U^ + Rj^T) - (U^ + SR^T)"^ + BZ2 


px / ^ \ • *' 2 ^ 

-n r . 9 

(uu - u") 


( 2 R 2 T‘^ + 2 R^T) 


+ BZ3 


•,;uu'*' (u"^ + 3R2T‘^) - u^(u'^ 


2 " 

+ Bziruu-" (n-^ + SBjT-*) + DD'"(U® + R^T) - 20^ (U'" ^P-zT*)] 



— 

BZ5 (UU'^ (2r?,2'^'^ + 2R^T) j- (s:) (x) + 

'.M- 

Hz f ■" 

(D-"^ + R^T*) - CU'" ( V*^ * SRiT*)] + BZsjTu^- UTJ-") 

( 2 B 2 T + 2B^T ■")'][ + BZ3 [u^ (U® -f SHjT) - TJV* (U^ ^ EjT) 


+ BZ4rn^ (U* + SRgT) + {n*% TI 1 T+) - 2CU* (C^^gT)^ 


EZsTu^ (2R2'^ 2R^T'^) I Kg (x) N^ (x) 


II 

X 22 ^ CZ3 (u^ + SRgT) - xnj-" (U^ + 

CZ4 r U^ (u^ + 3 R 2 T) + K^(u'' + ■" 

+ CZ5 r U^ ( 2 R 2 T +■ 2^2^"^) [ Ng Kg + 

. y 2 n 


^ 2 " 2 1 

X f CZ 3 u'^TJ (K'"^ + "■ ■" 

22 j 

CZ4 |^U'*’U (U’’" + 3R2T^) + U'^U (U^ + R2'^^ ” 

(U"^^ + R„T’^) + CZ5 j U'^U (2R2T'*' + 2 R 2 T) : Ng Kg 

^ "J (1.4 


( 1 . 43 ) 


where 



BZl = 

2 

(1 *Jif 

(3/.J - SAg + A 

BZ2 = 

3 

(1 *pf 

(A^ + ^2 - Ag) 

BZ3 = 

- 6A^/(1 + p) 



= 6 (2A^ - Ag)/ (1 +/?)^ 


BZ4 
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BZ5 = Zi\^/ (1 + pi 

CZl = - I (3/i^ - 3^2 ^ ^3^ 

CZ2 = - I (Aj + J'.j - -f. 3 ) 


CZ3 = - 3^.J^ 

CZ4 = I (2/i^ - ii'g) 

CZ5 = I ^2 
Z3 = 

Z4 =1 (2ii^ - A^) 


Z5 = ^ 

AZl = 2(3A. 


2^2 + ^3) 


A 


~\ 


AZ3 = 3(A^ + ^2 *■ ^' 3 ^ 

AZ4 = 6 ( 2A^ - Ag) I 


J- +y3 

1 +j-i 


A 


AZ5 = 3/i, 


/S 


2 ' 11 + A 
22 




1 + 


P 


and the d ^ -moment transfer equations for rigid 

X 

sphere molecular model. 




(t"^- T) - 



-■ k" nJ' e(sj) (T + t'^)J 

- N>; m/ \IV * 3V!j 


3/2 



as 


-M* N- m/ j 




[-1’ 


3/2 


(1.44) 


and 


W \ ^""^2 
d£ } 2 

f I 


3Ro (T""- T) - (U^ - U-" ) 


- i/ir-x n"n^ ©(sg) (2R2T ->'■ 


+ 3 /a 


n*n; JA I-;/ < n 2 ^ 'Is > 

- N^N* cq-g 0 (s„. , U)4--;>{4) 


2R^T^+ 2R2T 


2R^T + aRgT 


3/2 
3/2 

(1.45) 


where 


1/2 


1/2 


= (U - U'^)/ p2R^ (T + 

82 = (U - U"")/ j^2R2 (T + T'^)j 
©(s) = (s^ + 1 - — ^ ) exp (-s^) 


4 s 


. (2s3 * 3s - J- . -% ) I erf(s)/A 

4s 


M. 


= M^/ (m^ + Eig) 


M2 

^12 

isd^ 

ijOo 


VK “^2^ 


= (cTi + r 2)/2 




n. 


“a'^2 ( il * i; > 



T 





^ 2^-1 ^ 



®ud = (SRgT + 2R^T'") 


0{s,(M),ri) = j (Sfj - i)s^ + (2Ca)+ V]- |) 

+ (|-<«^) s^ + |^(C«J- |) + (4ix>-6+2i| )s 

s 


T. 

2 (V| - |-)s^ l^erf(s) 


Detailed study Of the problem for various values 
of concentration ratio, mass ratio and Mach Number is 
not reported in the literature. The main objective of 
the thesis is to discuss the nature of system 1.40, 1,41, 
1.42, 1.43, 1.44, 1,45 and 1.25 and to present their . 
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solutions for various sets of the values of (OC , 
to bring out - 

i) The individual effect of these parameters 

2 3 

ii) A consistent tJT^pd, if any, in the and c„ 

moment solutions and 

iii) The effect of molecular model on the shock 
structure. 
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In the present investigation, the shock structure 
is determined provided the molecular model and the two 
additional velocity moments of the Boltzmann equations 
are prescribed. 


The conservation of mass, momentum and energy yields 
the following governing equations; 



leading to the following expression for N^, Ng , 

, 4 * 


~ and 
u 


L 

T 


N- = 1 - N+ 


( 2 . 5 ) 


N; =: n; 


( 2 , 6 ). 
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ir 

u 


SL 

T 


1 

4 


(1 


5R^T 


U 




= 1 + 


\J‘ 


5£^T 


r -Tr-i 
1 1 4 2- 

L ■ ^ 


^1 \/^ ^ 2 


(2-7) 


( 2 . 8 ) 


For the remaining unkno^jns Nf and Nl we have to 

1 2 

use two extra moment (trf^nsfer) equations. As in these 
extra transfer equati onf , the c orresponding moment of the 
collision integral will not vanish, the molecular model 
will also play a part. To verify if the choice of the 
transfer equation and of the molecular model have a con- 
sistent trend in the final result s, we shall study three 
different cases; 


Case 1 , d - moment and Maxwell molecule 


-1 > Mi! 

d? ■] U 

/ I 


N] 


3R^T 


+ . 2 




i -■ 


^12 i 

I up. 




u _ u 


(OC - Np ) (l - -7 + 


'1 U 


U 


»2 - ® ^ TT^ f fl _ ^ 

- A 2 

[Sf] I ^2 ^ ® nJ - 2 2" n; 


2 , 


( 


- Agd - f Np n; (1 - 


U\2 


ft ^ > 

- U / j 

(2.9a) 
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j „+ . \ 3F,„T „ + 
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T 
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(|- - 1 ^ - (1 -I 


/ u *^ ^ F ) 


U 


’jj ' 


T+i ! f ® " §■ ^2 ^ ''i 1) + (1 - §- Np 

L r I — 


N, 




2 U 

(M. - 2/\„) 


(t ' (rr M+ - A - Z'^I 1 '4 ! 

(i - U ^ + (OC - N 2 ) 2 l+£ ^T ■* M 

u / — ‘ 


(N* + ffiN^ - 2 ^ nJ N* ) (1 - f)J 

+ / r - ) 


U‘"s2l 


(1 

/ 


U ,T+ \ 1T+ / J U \ 2 


- X22 (ce_ ~ Np n; (1 - 


u 


u 


(2.9b) 


Case 2, 6„ - monent and Max-sjell molecule : 

X 


^ i N - 


6R.T 3R.^T^ xi TT + 4 TT + 2 

( 1, _1_, _J:_) . 6 (|-) 

I 


U * 
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R^T ^ + 


u ‘ 


T 
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,+ p R.T 3R.T 
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- 4 - (1 + 1 -) 
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+2 R^T 
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2 T 
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T J 


(2/i. - Aj 


U 1 \ 

^ (1 + -g- ) 

U 


SlSfSll ^ 3 V t ; 

■*■ U 1 'TJ ^ ^ 1 
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Case 3# "* moment, rigid sphere gas 
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and 



where & and 0 are the same functions as defined in 
Chapter I. j. . 
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Each of the system (2,9), (2,10) and (2,ll) consists 
of a first rorder nonlinear, coupled differential equat- 
ions and is subject to the following boundary conditions. 

As 00 , nI" . 1 (2.12a) 

Y 1 



The solution is physically acceptable, if at the other 
endji.e.j^ — - oo , it satisfies 

Lim = 0, Lim = 0 (2,12b) 

OD 00 

It appears that we may start the integration from 
either end and verify that the conditions at the other 
end are fulfilled, H owever, in practice, the procedure 

is not as straight forward as it seems} one can start in— 

/ 

tegration from only that end which satisfies the criterion 
for being a ’saddle point, ' This point will be elaborated 
in the following. The detailed discussion will be given 
for the system (2,9), 

Theory: 

To analyse such differential systemswhich are auto- 
nomous (because there is no explicit df penyahce of inde- 
pendent variable tJn slope function) phase pjfene-methods are 
used. The singular points are located and the behayiour 
of the integral curves in the neighbourhood of the singu- 
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lar Doints are examined. The knowledge of the nature of 
singular point and the slope functions enables one to 
predict the qualitative features of the system. 

For example, let us consider the system. 


da 

6^ 






(2.13) 


This system is autonomous because there is no ^-depen- 
dence on the right hand side. Moreover Y = X = 0 is a 
singular point because both the slope functions vanishes 
for these values of J and Y. Looking for solution in 
the (X,. Y) plane, we write as 


51 - 

dX " 


a^l X + a^^Y 
^11 ^ ®'i2^ 


(2.14) 


It is possible to start the integration of (2,14) 
if we know the starting slope S, where 6 is the limit, 
as the singular point is reached, of the ratio y/X. 

The procedure to calculate S consists of first 
solving the following algebraic equation (in an unknown A 
- nothing to do with the mean free path) 

i 

^11 ” X ^12 =0 

®^22 ^22 “ X 


If the values of \ are real and of opposite signs, the 
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singular point is classified as ’saddle point'. If the 
values of a , are real and of same sign it is cl§.ssif led 
as ’node’. When the singular point is saddle there are 
only four integral curves passing through the point and 
the selection of one of then rj^hich _s relevant :t'0. the 
problem becomes possible, 

Gilbarg (1951) applied this method to study the 
shock wave structure in a single gas using the continuuta 
theory and proved that in the single gas case, as long 
as the gas satisfies a certain (fairly liberal) equation 
of state, the downstream singularity is saddle and the up- 
stream one is node in character. 


However, we shall show here below that such a gener- 
al result is not available for our (mixture) problem; 
the nature of the singular points depends on the set of 
values of (OC . /3,M), We start with the singular point at 
the downstream infinity. 


Downstream infinity ; 
Here ^ 

- T 

OC 


Y = - 


+ N" 


+ n; 


so that the linearised form of the equations (2,9 a-b) , 
valid in the neighbourhood of this point, are 
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Where 
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The roots of the characteristic equatio<j are and 
which satisfy 

^11 “ A ®12 

= 0 

^21 ^22 “ X 

The values of for various values of the set of 

parameter ( OC^ ^ ^ M) are given in the tahleX It is 




observed that in only a few cases I’s this end a saddle* 
For such set of values of (OC^ M) integration may be 
started from this end. For others we have to examine 
the nature of the upstream infinity. 


Upstream Infinity; 


Here we do not need to define new X and Y as: 
and N 2 themselves vanish at this singular point. How- 
ever, one cannot obtain directly the linear form (2,13) 
of the equations (2.9) because of the repeated occurrence 

“f" "4“ 

N + N / / 

of the expression 1 2 (which is of the form 0/0 


'I 


NJ +PnJ 


at the singular point) and its derivative. 
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1 ' 
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Consequently, when the system (2,9) is linearised with 
respect to N^j the elements of the coefficient matrix 

a will contain 5 , in addition to M. The final 


linearised form is 
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^12 ~ ^^2^ 
(1 


(1 - V*) - X 22 ^2 ® 


Eigen values of the coefficient m trix are given by 

^ 1 ’ ^2 “ i“ ^ + ‘^ (^) ) 1 ■! i ;- 

I J 


+ 4b(&) c(^) 


As the roots Aj[» A 2 characteristic equation 

of Jb involve A » one cannot determine their sign from 
the knowledge of OC , ^ ,1' only; one has to solve for S 
from other consideration and then come back to determine 

, f 

the sign of A s. ^7e proceed as follows; 

Assume (to be verified a poster iori )that the up- 
stream infinity is a saddle. Then •£ is determined as 
solution of the algebraic system 

aTty -A. U) 


a (5) -AgTAJ 

This algebraic system admits a finitely large number of 
solutions but the one satisfying the inequality 


-a (S 

"W 


< & < 


-e(t 


=if < S 


-a (A) . 


is the one we seek. 
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This leads to the following conclusion for system 
(2,9) expressed in the form; 

= f (<, Np 


dN 




dN" 




= g (N^, N^) 


Firstly if there are only two solutions of the system 

i) g (N^ , Ng ) = 0 

f (N^ , ) = 0 


Let us denote them hy Z^, 3^ 



ii) For a fixed in tie interval 

"+ 1 
0 ^ ^ |- 

s 

Where g (N^, r) = 0 and 
f (N^, s) = 0 

iii) In the region hounded hy the curves 
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0 


iv) The point is saddle in character 

and (v) There is only one curve having direction^at 

the saddle point which satisfiesthe inequality 


-f + 

^2 


s 




Z. 


Then there exists an integral curve joining and^which 
satisfies the diff erational equations and is unique. 

For such values of the set (OC , ^2 ,M) which make 

the upstream singularity a saddle point, integration is 

started from upstream end. There are values of the set 

(OC , , M) for which neither upstream nor dov/nstream 

happens to be a saddle and the corresponding solution 

cannot be obtained by the systematic procedure. Here 

one has to resort to a guessing procedure, i.e,, assume 

the values of at a station (say -^ = 0, = a where 

1 + 

0<a< — - I in particular N. was taken to be equal to 

^ in 

0.5 for numerical integration) and integrate^both direct 
ions, i.e. towards - ooand + oo . If the boundary condi- 
tions are not satisfied, change the assumed value of 
Ng and repeat the process till a reasonable integral 
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curve is obtained. We have obtained solutions for three 
set of values of M) by this guessing procedure. 

Discussion ; 

A shock wave passing through a homogeneous gas mix~ 
ture induces separation of species, i.e,, the concentration 
of species changes in the shock wave. In the present con- 
text we shall not use the term first and second species and 
concentration ratio will denote the ratio of the number 
density of heavy species to the number density of light 
species and relative concentration ratio will be defined as 
the ratio of concentration ratio to initial concentration 
ratio. In order to bring out 

a) The effect of initial concentration ratio o<(^, mass 

ratio and freestream Mach number M on the shock structure, 

2 3 

b) A consistent trend, if any, in the c and c moment 
solutions and 

c) the effect of molecular model on the shock structure - 
the computation of number density profiles was carried out 
covering a good range of concentration ratio, mass ratio 
and freestream Mach number. The results have been presentee 
graphically in Figs. 3-30 and are discussed below; 

Effect of concentration ratio ; 

The result of computation corresponding to|^= 10, 

M = 2 and ranging from ,022 to 2,5 has been reported in 
Figs, 3-6, 9 - 13, 15 and that of corresponding to|f= 10, 
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M = 3 and 2*0 to 3,0 has been reported in 

Figs, 19 - 21. 

It is seen in these figures that the relative density 
curves do not intersect below a certain value of concentrat- 
ion ratio for a given mass ratio and free— stream Itach number 
Thus 5 there seems to exist a critical concentration ratio 
for given pair of species depending upon ^ and M such that 
for any concentration ratio greater than the critical value, 
relative concentration ratio assumes a value greater than 
unity. This can also be expressed as ’the velocity of the 
heavy species remains greater than that of the light species 
throughout the shoclc wave as long aa the concentration ratio 
does not exceed the critical concentration ratio’. It can 
also be seen that for a given^and M the number density at 
which intersection occurs decreases with increase in the 
concentration ratio and the overall rate of compression of 
heavy species increases with concentration ratio. 

Such a trend of density profile is not found in any 
of the earlier numerical investigations. Bird (1968) has 
not made any computation in the range where this trend can 
be observed. Goldman and Sirovich (1969) have made computat 
ions in this range and have found the velocity of the heavy 
species higher than that of the light species throughout 
the shock wave along with the velocity undershoot of light 
species. Center (1967) claims to have done the experiment 
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in this range hut neither density profiles beyond cX. = *923 
have been presented nor any remark about the intersection of 
density profile has been made. The density profiles corres* 
ponding to(?^ = .923, ^ = 10 and M = 2.25 have been reproduc( 
in Fig, 5 along with Center’s experimental result. Inter- 
section of density profiles can be seen in these figures. 
The second source of qualitative information on this subjec 
is the asymptotic analysis carried out by Harris and Bienkow: 
(1971). Using asymptotic arguments they have shown that the: 
exists more th?^n one length scale of transition depending Oi 
the limiting values ofit^', , M and their products. Accord' 

’» r ' 

ing to their analysis each of oCj and M is of order 

one, the transition occurs on one scale, i*e*, light and 
heavy species are strongly coupled and the dissipative ffiec— 
hanisms resulting from self and cross— collisions are of ord 
unity* On the other hand, when both cX and are large 

in the limiting case on the inner scale, heavy-heavy self 
colllsional shock is formed and light species is constant 
to the zero order* The pattern resulting from Figs* 12—1 
seems to be In conformity with the view point of Harris and 
Bienkowskl on account of regular shift of the point of In- 
teraction of density profiles and regular increase In over- 
all compression rate of Jargon* Thus we conclude 

a) velocity of the heavy species is not always 
higher than that of the light species throughout the shock 


wave, 


5 : 


b) the velocity of species decreases monotonicallj 
to their corresponding Rankine-Fugoniot values 

and (c) the pattern emerging from rate of compress- 
ion of individual species and region where the species ar{ 
significantly compressed (e.g,, the region of significant 
compression of heavy species may be different from the re- 
gion of significant comnression of light species) suggests 
transition on multiple scales and formation of shock withi 
shock. 

Effect of Mach number ; 

The density profiles cor res bonding to M = 2 and M = 
3 have been given in figs, 9 and 18. The values of (J^and 
, the additional moment and the molecular model are lef 
unchanged and their details are given in these figures, I 
is seen that separation of species measured in terms of 
difference of their relative number density or, eqiiivalent 
the difference of reciprocal of their bulk velocities in- 
creases with increase of freestream Mach number. The figu 
also suggest that heavy species starts getting deccelerate 
anly after the light species has been compressed signifi- 
cantly. This trend persists in c moment solution also 
as can be seen from Fig, 27, Examination of Figs, 13 and 
19 shows that critical concentration ratio increases with 
increase in Mach number. 
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• Effect of mass ratio : 

The density profiles corresponding to ^ = 10 and 
= 2 have been given in Figs. 11 and 17, The other para- 
meters are common to these sets and are given in the figures 
As expected, more dissimilarity in mass creates more separat 
ion betnveen the species, /'nother aspect of the effect of 
mass ratio can be seen in figures 12 and 16. These figures 
imply that as the mass ratio increases, while other para- 
meters are held constant^critical concentration ratio de- 
creases. It is also observed that the effect of mass ratio 

that 

is more pronounced than/of concentration ratio on the se- 
paration of species as can be seen from Figs. 24 and 27 in 
which the effect of decreasing the concentration ratio by a 
factor of 10 gets obviated by decreasing the mass ratio by 
a factor of 2. 

Effect of additional moment on solution ; 

Though the overall qualitative asoects of the solut- 
ion do not change with the choice of additional moment, 

a sharp quantitative difference can be seen from Pigs, 6 
3 

and 22. C moment solution gives lower critical concen- 
tration ratio and it gives poor agreement with available 
experimental results. 
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Effect of Eiolecular model on solution ; 

Figs, 29 and 6 can Ise compared for this purpose. In 
this case, compression of both the species for rigid sphere 
gas takes place at a lo;i?er rate, thereby, thickening the 
shock. Separation of species is greater in mole- 

cule model. For rigid sphere gas, increase of Mach number 
increases separation bet;?een the species and delays the 
compression of heavy species in the beginning also, as can 
be seen in Fig. 30, It may be aporopriate to point out that 
with Mott-Sraith ansatz Maxwell molecule model yields better 
agreement with available experimental results. 

Effect of force constants ; 

The situation to discuss this would not have arisen 
had the traditional values of force constant deter- 

mined from data of diffusion coefficient (based on near- 
equilibrium experiment) and Chapmanr-Enskog theory (whose 
validity is questionable in strong-non— equilibrium situat- 
ion) - yielded satisfactory results. Phenomenon exhibiting 
diffusive effects being more complicated than the pheno- 
menon exhibiting viscous effects, numerical experiments 
were carried out to see the effect of density pro- 
files. It so happened that the change of from 1.385 to 

1,0 improves agreement with all the three available experi- 
ments. The effect of change of density profile can 
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be seen in Figs, 6 p.nd 7 which also suggests that increase 
of X ^2 increases the overall compression rate of heavy 
species. 

Sensitivity of numerical scheme : 

Though the numerical difficulties in obtaining the 
solution arc enormous, it is Interesting to note that the 
fin.al result is almost insensitive to higher order inte- 
gration scheme. First and 8th order integration schemes 
with variable step size were used to compare the density 
profiles and they have beenreported in Figs, 13, 14 for 

2.0, ^ - 10,0 and M = 2 and in the Figs. 7 and 8 for 
C< = 0,1, = 10 and M = 2, It seems likely that the 

effect of higher order Integration is counterbalanced by 
the round-off error resulting from lengthy and ill-condit- 
ioned numerical operations. 

Comparison with experiments : 

1 Comparison with experiments of Center (1967) and hiuntz 

(1972) have been given in Figs. 3-5 which also contain the 
details. The agreement with Harnett and Muntz’s experiment 
which corresponds to M = 1.58, ^ = 10, c< = 0,111 is very 
good, as can be seen in Fig. 3, Figs, 4 and 5 present the 
comparison with Center’s experimant which correspond to 
M = 2,07, (3 = 10 and c< = 0.022 and M = 2,24, (3= 10 and 
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CX, = 0.923. It is seen that in the forner ease separation 
of snocies is slightly overpredicted over a small interval 
and in the latter case the separation of species is under- 
predicted over a large oart of the shock wave. Inspite of 
the qurntitative agreement with experiment not being very 
»lose, as seen in Fig, 5, it is only the present computation 
which agrees with experimental results as regards qualitative 
aspect of intersection of density curves. 

Conclusion ; 

Thus the numerical experiments carried out lead to the 

following conclusions: 

2 

i) The Eioiaent and the model of Maxwell molecule 

X 

give the most satisfactory results, 

ii) Separation of species depends strongly on the model 
of molecules and the additional moment chosen, and 

lii) Intersection of density curves of heavy and light 
species defends on the concentration ratio, mass ratio and 
free stream Mach number. Patterns resulting from the variat- 
ion of these parameters suggest that the f ormation of heavy- 
heavy self GOllisional shock on the inner scale and transit- 
ion of light species on the outer scale is possible for a 
large number of cases characterised by the limiting values 
of these parameters. This trend is in fair agreement with 
the asymptotic analysis of Harris and Bienfcowski (I97l) and 
the experiment of Center (1967) (Fig. 5) and is not found in 
earlier numerical investigations, 
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The arbitrariness in the Mott-Smith method (namely, 
the choice of the extra moment) and the criterion first 
set up by Oberai (1967) to remove the arbitrariness was 
discussed in the introduction. Tfe shall extend his work 
to the mixture problem. The problem will be formulated 
for arbitrary Mach number, keeping the extended Mott- 
Smith aansatz Intact. 


Thus, the postulated distribution functions are 
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leaving only the two functions (x) end to 

be determined. 


The genesis of the method is th^t the unknown 
functions and may be so determined that the equat 
ions 
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’^ith the form of given in (3.1) the integral and 
^2 evaluated for either (Maxwellian and hard 

sphere )molecular model. Sowever I„ involves the inte- 

O 
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gration of which, in general, is a gigantic task. 


However, as pointed out hy Oberai (1967), in the 
limit: upstream M->oo , the supersonic modes in fj^ 
approach 'Dirac's delta function and Ig can be evaluated 
analytically for the hard sphere model. 
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To write the above mentioned result, we note the follow- 
ing integral. 
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Introducing a new function 
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the collision Integr.n can be written as, 
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l¥e want to determine tJ- and n' subject to mlnlniz-tlon 


Of the tinctional 
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UT = P'l (1 - n: (oc- n:) u(aB)] 

J s/a-*- 3/2 2 I J 

(h) [H%I 


+ 0(1!^) 


Finally it amounts to minimising tlie integral 
dN“ 




-t 2 


a«i r -I 

2 ^ £n“(1 - N") MO'..) + N“(0C- Np M(OCp)J 

(1 - Np M(i^.) + (oc - Np M(a/5> 

2 ^ In” (1 - nPm(. 0 + Ng (® - Np M(A^) J + 

%•* rt 


J^N“ (1 - Np M(/ ) + nJ (oc - Np M(aI^)J 


dx (3.23) 


This yields the following equations, 
dN’ 
dx L dx 

dN^r r t A 

(1 - 2Np M(/i) + Npl - Np M (A) + 


+ I^Npl - Np M(A) + NpOC- Np M(Ap)j| 


dx 


[ 


5n, 


(oc - Np M(/i p) + N” (oc - Np m'(^. p) + 

(3 N 

( 1 - Nj) m'(/i) ^ **2 <®" 

'U p) ^ (1 - Nj) M(Ji) + (CC- Np 


M 

M(Ap)J 


^ Np L 



(1- 2Np ?i(. ) + (l-Np 0^.) ^ (cc_ 


• i.. v-'-; — ^ \,uj - N2 )m(A|3) + 


N"(fE ~ Np ?” (;. p) p> 


No(® - N^) I"' ) 


fVt 

CM 


n 


5nu 


+ j N^Cl - n7) M(A) 4 


L. 


No (1 - nT) I.<a U) 


^n' 






1 

(3.24a) 


and 


dx 


dN. 


~35E '*' ^^2 ~ ^ 1 ^ r. (is) + 1^2 (OD - N 2 ) M(Ap) 


d N o I « •— A 

(I - N“) K(A) 4 NgCl - N^) i: (A) ^ 4 


dx 


(0C-'2Np L(A/.-- ) + N" (OC - Np M 


dNp 


dN< 


dx 


n“( 1 - Np hpA) 


M_ 

^n: 


- N 


“ M(Ap) 4 NpOC- Np 


M- 

I 

K(AiG ) 


^>A 


^N" 


Np 1 - Np M(A) 4 NpOC - Np 


(1 - Np N(/) 4 N 2 (1 - Np M'’(A) 


dA 

diC 


11 




c)n: 


(OC- 2Fp H(A(i) 4 NpOC - Np M*(Ap) 


iq (1 - Np H(A) 4 N; (OC - Np M(Ap 


n7(i - nT) ?r’(A) + 


Da 


■1^^ - 
Np(A|3) 


Bn7 


N-(a- Np uHAfi)f — 


(3.24b) 
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TItf'sr I’quations have to be solved subject to the 
follo'A'iJifT hmindary conditions: 


Lim 

X"»00 




Lim n't = I 

X-*. - CO 



3.24(c) 


The system (3.24) is far more complicated than that 
encountere*} in Chapter II, •That u'e have gained in terms 
of freedtJn from arbitrariness is more than offset by 
the almost unnanapeable complexity of the (differential) 
system to Ir solved. 

For a for typiC'^l cases the numerical so?.utions 
obtained by three-fold guessing (.i.e. at x = 0, one 
chooses a reasonable velue of n 7 - > sa;^ 0,5>and guesses 
the corresponding values of II,,, that lead to 

correct behaviony at + oc^re given in Fig, 

In conclusion, it may be mentioned that for the 
bimodal function (i.e. for the single gas problem) 
Narasimha and Desbpande (1969) have succeeded in evaluat- 
ing, in closed form, the collision integral J, for the 
hard sphere molecular model. Tiowever, the related inte- 
grals appearing as co-efficients in Euler’s equation do.,3 
not depend on position and this simplicity is lost in the 
corresponding problem for the gaseous mixture and the 
problem remains untractable. 
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APFSNDir. 


*J^wo of Rtftssos £iti(3 ^2 iiioving Tiritli 

volocitiGs Cji^ fl.nd col lid© ftnd sir© tlicir finsl 
velocities. The initial relative velocity is devoted by 
g and J is the unit vector along the direction joining 
the position of the particle and the point of inter- 
section of the asymptotes to the trajectory, describing 
the relative motion of m2 (Pig. 


Then (cf. Ref.( g), Sec.3.43^ chapman and 

®1 ~ £1 = 2(mV(m + m ) (g, j) j 
~ - 2(m./(m. + m ) (g, j) j 

MM JL X A 


Cowling 

(1964)) 


Again (see Pig. 1(h)) 

1, ^ = Cos © Cos]4^ + Sin © Sini/i^' Cos 0 

therefore 

= (c^ - £l)* i 

- 2 g Cos © (Cos © CosiJw + 

+ Eig; ° i 

Sin 9 Sinf: Cos 0 

f “1 i • 

Uo - Uo = - 2 ■ g Cos © (cos © Cosy^ + 

A 2 KS^ + / 

Sin 0 Sin^' Cos © 

Also cos^ = (u2 - u^)/g 

and the angle of deflection '])c =1T- 2 © 



Calcul^ion of 


J 0 (c) J (f, f) for Maxwell molecule require 
evaluation of 


00 

f r - V ^ 

(l - Cos denoted by Aj 

J ^ 

o 

which has been taken from GhfKp^n e^d^iGhw^lng? (1964) and 
Mott-Smlth (1954). 


TABLE - I 

DOWNSTREAM SINGULAR POINT EIGENVALUES X i 



- 86.9460 - 85.9564 



































TAB3IE II 

UPSTREAM SINGULAR POINT EIGENVALUES X i 












TiBLE.II (contd) 

UPSTREAM SINGUIjAR POINT EIGENVALUES X -j end 
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